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Background/Aims 

Small intestine motility is governed by an electrical slow wave activity, and abnormal slow wave events have been associated 
with intestinal dysmotility. High-resolution (HR) techniques are necessary to analyze slow wave propagation, but progress has 
been limited by few available electrode options and laborious manual analysis. This study presents novel methods for in vivo 
HR mapping of small intestine slow wave activity. 

Methods 

Recordings were obtained from along the porcine small intestine using flexible printed circuit board arrays (256 electrodes; 4 
mm spacing). Filtering options were compared, and analysis was automated through adaptations of the falling-edge varia- 
ble-threshold (FEVT) algorithm and graphical visualization tools. 

Results 

A Savitzky-Golay filter was chosen with polynomial-order 9 and window size 1.7 seconds, which maintained 94% of slow 
wave amplitude, 57% of gradient and achieved a noise correction ratio of 0.083. Optimized FEVT parameters achieved 87% 
sensitivity and 90% positive-predictive value. Automated activation mapping and animation successfully revealed slow wave 
propagation patterns, and frequency, velocity, and amplitude were calculated and compared at 5 locations along the intestine 
(16.4 ± 0.3 cpm, 13.4 ± 1.7 mm/sec, and 43 ± 6 |iV, respectively, in the proximal jejunum). 

Conclusions 

The methods developed and validated here will greatly assist small intestine HR mapping, and will enable experimental and 
translational work to evaluate small intestine motility in health and disease. 
(J Neurogastroenterol Motil 2013;19:179-191) 
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Introduction 

A number of coordinated regulating factors contribute to 
governing the motility of the small intestine, including neural, 
hormonal and bioelectrical systems. One of these factors is an un- 
derlying rhythmic electrical slow wave event that is generated and 
propagated by a network of interstitial cells of Cajal (ICC). These 
slow waves coordinate the timing, frequency, and direction of in- 
testinal contractions.^ However, the mechanisms contributing to 
normal and abnormal small intestine slow wave pacesetting and 
propagation remain unclear and controversial, largely due to limi- 
tations of previous recording and analysis techniques. 

Previous studies have measured slow wave activity in the 
small intestine, but most have been restricted to sparse recordings 
where a few electrodes are spaced along the intestine. ^'"^ Such 
methods provide information about slow wave morphology and 
frequency, but they are unable to provide spatial information 
about slow wave propagation dynamics, such as pacemaker activ- 
ity and wavefront interactions.^ High-resolution (HR) techni- 
ques have more recently been introduced, whereby recordings are 
obtained from dense arrays of many electrodes, allowing activa- 
tion sequences to be recorded in fme spatiotemporal detail.^ ^ 
These HR mapping techniques have now been applied to eval- 
uate gastric slow wave patterns in health and states of dysrhyth- 
mia, including in human motility disorders. ^'^ Clinical interest in 
evaluating intestinal slow wave patterns has been fuelled by asso- 
ciations between a range of motility disorders and aberrant slow 
wave activity or ICC function, such as diabetic dysmotility,^ irri- 
table bowel syndrome, and post-surgical dysmotility.^^ Howev- 
er, there remains a critical lack of clinical intestinal slow wave 
studies to date. 

While an automated pipeline of sophisticated algorithms now 
exists to support gastric mapping, ^^"^^ there is no equivalent sys- 
tem for analyzing the vast data sets that result from HR small in- 
testine mapping. This data throughput problem has been a major 
impediment to experimental and translational progress. In addi- 
tion, small intestine HR mapping to date has been limited to flat, 
rigid electrode platforms. ^'^ While these arrays achieve a high sig- 
nal quality, they cannot conform to the intestinal curvature, so 
they can only be applied to a portion of the circumference. ^'^^ 

The aims of this study were therefore: (1) to develop reliable 
HR recording methods for in vivo intestinal electrical activity us- 
ing flexible printed circuit board (PCB) arrays and (2) to im- 
prove the efficiency of intestinal slow wave data analysis through 



the development of automated algorithms. The practical utility of 
these advances was validated in a porcine model. This work 
builds on that initially presented in abstract form at the IEEE 
EMBS conference. In particular, the effects of 4 different fil- 
tering methods were quantified, intestine slow wave velocity and 
amplitude maps are presented, and a more complete set of ex- 
perimental data was included in the optimization of the analysis 
algorithm parameters. 

Materials and Methods 

Electrode Design and Application 

Flexible PCB arrays were designed for intestinal recordings, 
based on those validated for gastric mapping.^ These arrays con- 
sisted of gold contacts with copper connectors, embedded in a 
polyimide base. The recording head of each PCB consisted of 32 
electrodes (0.3 mm diameter) arranged in a 16 X 2 array, with an 
inter-electrode spacing of 4 mm. The 4 mm spacing was calcu- 
lated to be sufficient for recording small intestine slow waves 
based on the spatial variant of the Nyquist theorem. Eight 
PCBs were used simultaneously in each recording, comprising 
256 electrodes in a 32 X 8 array, covering 124 mm of intestinal 
length and approximately 35 cm^ (Fig- lA). The PCB arrays 
were contained in gauze-padded silicone cradles that were sized 
to match the circumference of the porcine jejunum (Figs. IB and 

IC) . The cradled arrays were placed such that the electrodes con- 
formed to the intestinal curvature, covering most of the circum- 
ference, while exerting only gentle pressure on the serosa (Fig. 

ID) . Small windows were made in the mesenteric attachments at 
the edges of the cradles, through which gauze was tied over the 
top of the cradles to gently hold the electrodes in contact. 

Experimental Recordings 

Ethical approval for this work was obtained from the Univer- 
sity of Auckland Animal Ethics Committee (R698). All experi- 
ments were performed in vivo on white, cross-breed, weaner pigs 
of either sex (n = 5; 34.9 ± 0.7 kg). Animal preparation and an- 
esthesia were as previously described in another study evaluating 
gastric slow wave activity. In brief, the animals were fasted over- 
night and then subjected to general anesthesia, which was in- 
duced with Zoletil (Tiletamine HCl 50 mg/mL and Zolazepam 
HCl 50 mg/mL) and maintained with isoflurane (2.5-5% with an 
oxygen flow of 400 mL within a closed circuit anesthetic system). 
Slow wave activity has been successfully recorded under anes- 
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Figure 1. In vivo application of electrodes. (A) Eight printed circuit board (PCB) electrode arrays were tessellated to create a large rectangular 
electrode array with 256 electrodes (arranged in 8 by 32 configuration with 4 mm inter-electrode spacing). (B) Exteriorized small bowel segment. (C) 
Electrode array layered in gauze padded siHcon cradles. (D) Small bowel segment held in contact with the electrode array by the siHcon cradle. 



thetic in many previous studies in several species, ' ' and anes- 
thesia of this type is thought not to have a major effect on slow 
wave activity. 

Following midline laparotomy, the flexible PCBs were ap- 
plied to the serosa at intervals along the small intestine. The por- 
cine small intestine was approximately 15.5 m in length, and the 
first recording was obtained from the duodenum, immediately 
distal to the pylorus. Two jejunal recordings were taken, the first 
immediately distal to the ligament of Treitz (approximately 0.5 m 
distal to the pylorus in pigs; proximal jejunum) and the second at 
4 m distal to the ligament of Treitz (mid-jejunum). Two ileal re- 
cordings were then obtained, the first was at 9 meters distal to the 
ligament of Treitz (proximal ileum) and the second was at 1 0 cm 
proximal to the ileo-cecal junction (distal ileum). Visceral hand- 
ing was kept to a minimum. The intestine was placed back into 
the abdomen during recordings, and the wound edges were ap- 
proximated and packed with warm wet gauze. Unipolar record- 



ings were acquired at a sampling frequency of 512 Hz using an 
ActiveTwo System (Biosemi, Amsterdam, The Netherlands), 
with reference electrodes on the hind leg. 

Filtering 

Raw data was down-sampled to 30 Hz for computational ef- 
ficiency, and baseline wander was removed using a Gaussian 
moving median filter with 20 second moving window. The fre- 
quency of interest was determined to be in the 0-120 cpm (0-2 
Hz) range by analyzing the frequency domain via a Fourier 
transform of the baseline-removed signal (Figs. 2 A and 2B). The 
performances of 4 digital filters were investigated, with parame- 
ters set to accommodate the frequency of interest and as pre- 
viously set for gastrointestinal extracellular studies: Savitzky- 
Golay filter (polynomial-order 9, window size 1.7 seconds), 
Butterworth band-pass filter (5-60 cpm), Butterworth low-pass 
filter (120 cpm), and continuous wavelet transform (CWT) 
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Figure 2. Filtering of intestinal slow waves. (A) The frequency 
domain of the raw unfiltered signal shows high levels of low- frequency 
noise. (B) The frequency domain of the baseUne-removed signal shows 
that the moving median filter effectively removes the low-frequency 
noise and the frequency of interest lies in the range of 0-120 cpm (0-2 
Hz). (C) A representative raw slow wave signal from one channel of the 
printed circuit board electrode (a) is compared to the filtered form using 
the 4 selected filters: Savitzky-Golay (b); Butterworth 5-60 cpm (c); 
Butterworth 120 cpm (d); and CWT 5-60 cpm (e). The Savitzky- 
Golay filter maintained the best signal quaHty (see also Fig. 4). 



band-pass filter (5-60 cpm). 

Slow wave activation times (ATs) were manually marked in 
30 unfiltered channels from 3 separate recordings (n = 3 pigs). 
Manual marking was performed by visually determining the 
point of steepest negative down-slope on the slow wave, which 
corresponds with the arrival of the wavefront under the recording 
electrode. Using these manually-marked ATs, 3 metrics were 
computed for evaluating the filters: noise correction ratio (NCR), 
slow wave amplitude and slow wave gradient. NCR estimated the 
level of noise removed in sections of the iso-electric regions of the 
leads, computed as, 

llnunfiit - nfiitib 



NCR = 



|nunfilt||2 



where nunfiit is the noise of the iso-electric regions of the unfiltered 
signal, nfiit is the noise of the iso-electric regions of the filtered sig- 
nal, both after baseline removal, and H-.-l^is the Euclidean norm. 



Amplitude was computed by placing a window of 10 sample 
points (i.e., one-third of a second) on either side of the AT, and 
taking the difference of the maximum and minimum values. The 
gradient of the slow wave down-stroke was computed using a 3 
point central difference. NCR, amplitude, and gradient were 
compared between the filtered and unfiltered signals, to de- 
termine the optimal filter for removing low and high frequency 
noise while maintaining signal integrity. 

Signal Analysis 

Data analysis was performed in the Gastrointestinal Electri- 
cal Mapping Suite (GEMS) vl.5.^^ An automated marking al- 
gorithm (falHng-edge variable-threshold [FEVT] ) was adapted 
and tuned for marking the small intestine ATs.^^ In brief, FEVT 
calculates a specified signal transform and convolves the trans- 
formed signal and a falling-edge kernel. It then utilizes a 
time-varying threshold for detection of times when the convolved 
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signals exceed the threshold, thereby identifying ATs at the point 
of steepest negative descent. The FEVT algorithm includes sev- 
eral tunable parameters, including a threshold multiplier (r|), re- 
fractory period (Tr), smoothing kernel width (p), and running 
median window half- width (Thw). The Tj defines the threshold 
which the signal must surpass in order to be considered as a slow 
wave event, Tr defines the minimum period between FEVT 
marks, p defines the time width of the moving average filter ap- 
plied to the detection signal transform, and Thw defines the time 
window for which the time-varying threshold is computed/^ 
Additionally, the signal transform used for the FEVT algorithm 
was also optimized between 4 different methods: negative de- 
rivative, amplitude-sensitive differentiator, non-linear energy op- 



Table 1. Falling Edge Variable Threshold Algorithm Parameter 
Tuning 



Parameter 


Values tested 


Optimal value 


Signal transform 


ND, ASD, NEO, DE04 


NEO 


T] (unit-less) 


2,3,4 


2 


Tr (sec) 


0.5, 1, 1.5,2,2.5 


2 


p (sec) 


0, 1,2 


0 


Thw (sec) 


10, 15,20 


10 



A range of values of the FEVT signal transform were systematically combined 
and tested as listed in column 2. The parameter values that resulted in the 
combined best performance across all 3 data sets are shown in the third column 
and yielded performance metrics as presented in Table 2. 
ND, negative derivative; ASD, amplitude- sensitive differentiator; NEO, 
non-linear energy operator; DE04, fourth-order differential energy operator; 
FEVT, falling edge variable threshold. 
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Figure 3. Isochronal activation time (a), velocity field (b) and ampHtude (c) maps of small intestine slow wave activity. The electrode array was wrapped 
around the circumference of the jejunum such that the top and bottom of each map correspond to nearly-adjacent tissue at the mesenteric border. Three 
propagation examples are presented: (A) Consistent antegrade propagation (isochronal interval 0.5 seconds; see also Video 1). (B) Consistent 
retrograde propagation (isochronal spacing =1.0 second; see also Video 2). (C) CoUiding slow wave fronts. Slow wave activity originates both oral and 
aboral to the array and coUides in the middle (dashed Hne); (isochronal spacing =0.25 seconds; see also Video 3). The activation maps (a) show a single 
wave front, with each color band indicating the area of slow wave propagation per unit time, progressing from red (early) to blue (late). The velocity 
field maps (b) show slow wave propagation direction at each electrode with an arrow, and the magnitude of velocity is displayed as a color gradient, 
according to the scale shown. The amphtude maps (c) display the slow wave ampHtude as a color gradient according to the scale shown. Each black dot 
represents an electrode, and white dots outHned in red represent electrodes where activity was interpolated. 
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erator, and fourth-order differential energy operator. The neg- 
ative derivative method computes the first-derivative of the signal 
and considers only those values which correspond to a negative 
deflection in the signal. The latter 3 of these signal transforms can 
essentially be categorized as the first, second and fourth order dif- 
ferential energy operators. The transformed signal accentuates 
high-amplitude, high-frequency components of the original 
recordings. In general, the nth order energy transform is propor- 
tional to AnCOn, where A and 0) are the instantaneous amplitude 
and frequency of the original signal, respectively.^^ 

FEVT parameter values and signal transforms were in- 
crementally combined to analyze 3 representative data segments 
from the porcine jejunum (n = 3 pigs), with values as shown in 
Table 1 . FEVT results were compared to the manually-marked 
data sets to determine true positive (TP), false positive (FP) and 
false negative (FN) AT results. Hand-marked data were in- 
dependently verified by 3 experienced markers to ensure accu- 
racy, and only marks that all 3 markers agreed upon were used. A 
TP result corresponded to the FEVT mark placed within one 
second of a hand-marked event, a FP result corresponded to an 
erroneously placed FEVT mark, and a FN result corresponded 
to a falsely omitted FEVT mark. Three performance metrics 
were calculated based on the TP, FP and FN results, as follows: 



TP 



Sensitivity = ^ 

Positive Predictive Value (PPV) = - 
Aroc = Sensitivity X PPV 



TP 



(TP + FP) 




Video 2. Slow wave propagation in the retrograde direction. Each pixel 
represents an electrode and lights up when the slow wave passes that 
electrode. Slow wave activity propagates from aborad (right) to orad 
(left), and sequential slow waves are colored differently for easy 
differentiation (see also Fig. 3B) (Video clip is available at http://www. 
j nmj ournal. org/j ournal/j ournal_view.html .?year =2013 &vol =01 9&num 
= 02&page=179). 




Video 1. Slow wave propagation in the antegrade direction. Each pixel 
represents an electrode and lights up when the slow wave passes that 
electrode. Slow wave activity propagates from orad (left) to aborad 
(right), and sequential slow waves are colored differently for easy 
differentiation, (see also Fig. 3 A) (Video clip is available at http://www. 
j nmj ournal . org/j ournal/j ournal_view. html .^year =2013 &vol =019 &num 
= 02&page=179). 




Video 3. Colliding slow wave activity. Each pixel represents an electrode 
and lights up when a slow wave passes that electrode. This recording 
shows 2 slow wave fronts, one travelling antegrade and the other 
retrograde, colliding and terminating in the middle of the animation 
frame (see also Fig. 3C). Separate slow wave fronts are colored 
differently for easy differentiation. (Video clip is available at http://www. 
j nmj ournal . org/j ournal/j ournal_vie w. html .?year = 2 0 1 3 & vol = 0 1 9 & 
num = 02&page= 179). 
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Visualization of Slow Wave Activity 

Experimental data analysis was performed to show ex- 
perimental validation and proof of utility on 24 total recordings, 
encompassing greater than 11,000 individual slow wave ATs, 
from 5 locations along the intestine (n = 6 pigs): proximal duo- 
denum, proximal jejunum, mid-jejunum, proximal ileum and ter- 
minal ileum. All automatically detected ATs were manually re- 
viewed to ensure accuracy. An automated grouping algorithm 
(^REGROUPS') was applied to cluster individual slow waves in- 
to coherent wavefronts,^"^ and manually reviewed. Automated iso- 
chronal mapping of ATs was then performed to display and 
quantify the propagation of successive wavefronts in isochronal 
time maps,^^ as shown in Figure 3. Animations were also gen- 
erated to visualize slow wave propagation as flashing pixels, rep- 
resenting ATs occurring at each electrode, as shown in Videos 
1-3. 



Analysis of Slow Wave Characteristics 

Slow wave frequency, amplitude and velocity were calculated 
on the same 24 recordings described above. Frequency was cal- 
culated by determining the time between subsequent slow waves 
and is reported in cycles per minute (cpm). The amplitude of 
each slow wave was determined by detecting the peak and trough 
using a zero-crossing approach within a 1 .5 seconds window cen- 
tered at the AT, and the velocity was calculated by using a vali- 
dated smoothed finite difference approach. ^"^'^^ Frequency, am- 
plitude and velocity were compared between antegrade and retro- 
grade propagating slow waves. Automated velocity field and am- 
plitude maps were generated to visualize these characteristics in 
spatial detail (Fig. 3). Frequency, amplitude, and velocity are re- 
ported as mean values with standard deviation (SD) or standard 
error of the mean (SEM). 
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Figure 4. Signal comparison of 4 filters previously used for gastroin- 
testinal extracellular recordings. Slow wave amplitude (A), slow wave 
gradient (B) and noise correction ratio (NCR) (C) were compared. The 
Savitzky-Golay (SG) filter with polynomial-order 9 and window size 1 .7 
seconds performed very similar to the low-pass Butterworth filter of 120 
cpm, out-performing the band-pass Butterworth filter of 5-60 cpm and 
the continuous wavelet transform filter of 5-60 cpm in regards to 
maintaining slow wave amplitude and gradient. The SG filter proved to 
be the best filter for small intestine slow waves. The Butterworth filter of 
5-60 cpm achieved the best NCR, and may prove useful for recordings 
that have large amounts of noise interference, but it must be noted that 
the slow wave amplitude and gradient are appreciably decreased with this 
filter. Butter, Butterworth; CWT, continuous wavelet transform. 
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Results 

Experimental Recording 

The flexible PCBs were successfljUy utilized to record in vivo 
small intestine slow wave activity, yielding readily identifiable 
slow wave events, as shown in Figure 2C. They covered a length 
2.5 times greater than previous studies (124 versus 46 mm).^'^^ 
The silicone cradles efiectively conformed the flexible PCBs 
around the intestinal circumference and held the electrodes in 
contact with the serosa, yielding slow wave recordings in a mean 
of 75 ± 9% of the channels. 

Signal Processing 

The Savitzky-Golay filter with polynomial-order 9 and win- 
dow size 1 .7 seconds best maintained the slow wave amplitude, as 
shown in Figures 2 and 4. This filter maintained 94% of the am- 
plitude of the raw signal, compared to 89% with the Butterworth 
120 cpm filter, 66% with the Butterworth 5-60 cpm filter and 
71% with the CWT filter. The down-slope was more sig- 
nificandy affected by all filtering methods because these filters are 
designed to remove high-frequency content and smooth the sig- 
nal, but the Savitzky-Golay filter once again performed well, al- 
though slightly below the Butterworth 120 cpm filter. The 
Savitzky-Golay filter maintained a gradient of 57% that of the 
raw signal, compared to 61% with the Butterworth 120 cpm fil- 
ter, 33% with the Butterworth 5-60 cpm filter, and 31% with the 
CWT filter. Although the Butterworth filter of 5-60 cpm per- 
formed poorly in maintaining the slow wave amplitude and 
down-slope, it vasfly outperformed the other filters in noise cor- 



Table 2. Performance Metrics of Tuned Falling Edge Variable 
Threshold Algorithm 



Performance metric 


Value after tuning of FEVT 
parameters (mean ± SD) 


Aroc 


0.79 ± 0.08 


Sensitivity 


0.87 ± 0.04 


PPV 


0.90 ± 0.07 


TP rate 


0.87 ± 0.04 


FN rate 


0.13 ± 0.04 


FP rate 


0.10 ± 0.07 



Average falling edge variable threshold (FEVT) performance over the 3 
validation data sets using the optimal FEVT parameters presented in Table 1 . 
Results exclude channels lacking slow wave signal content. 
Aroc, area under the receiver operating characteristic curve; PPV, positive 
predictive value; TP, true positive; FN, false negative; FP, false positive. 



rection, with an NCR of 0.78, followed by the CWT filter at 
0.29, the Savitzky-Golay filter at 0.083, and the Butterworth 120 
cpm filter at 0.078. 

Automated Analysis Algorithm Optimization 

The FEVT algorithm parameters were optimized and vali- 
dated against a total of 7,525 manually-marked and validated 
slow wave ATs over 3 data sets (n = 3 pigs). Since sensitivity and 
positive predictive value (PPV) are inversely correlated in prac- 
tice, area under the receiver operating characteristic curve (Aroc) 
served as the best overall measure of FEVT performance. A com- 
parison of the Aroc values for each combination of FEVT parame- 
ter values across all 3 data sets yielded the optimal set of parame- 
ter values of: non-linear energy operator signal transform, T] = 3, 
Tr = 2 seconds, p = 0 seconds and Thw =10 seconds (summar- 
ized in Table 1). 

Channels that did not contain any slow wave signals, likely 
due to poor contact between the electrode and intestinal serosa, 
were removed from analysis. The performance of the FEVT al- 
gorithm with the above optimal parameter settings demonstrated 
an average PPV of 0.90 and Sensitivity of 0.87, which yielded an 
average Aroc value of 0.79. Additional performance metrics are 
summarized in Table 2 and illustrated in Figure 5. 

Visualization of Slow Wave Propagation 

Activation maps and animations were constructed to visual- 
ize the propagation patterns of intestinal slow waves and show 
proof of utility of the developed methods. The results of the 





I 60 sec 1 

Figure 5. Five electrograms of slow wave activity from the porcine 
jejunum with automated slow wave identification by the falling edge 
variable threshold (FEVT) algorithm represented as colored marks. 
Red circles represent a true positive mark, green triangles represent a 
false positive mark, and blue squares represent a false negative mark. 
Amplitude is automatically scaled for each individual channel, but 
averaged 32 ± 6 |LlV in the mid-jejunum. The 5 electrograms shown here 
are from adjacent electrodes, such that slow wave propagation across the 
channels can be seen, as represented by the arrow. After tuning of the 
FEVT parameters across 3 data sets from 3 pigs, the algorithm achieved 
90% positive-predictive value and 87% sensitivity. 
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Figure 6. Average slow wave characteristics from five locations along the small intestine: proximal duodenum, proximal jejunum, mid-jejunum, 
proximal ileum and terminal ileum. Frequency (A), amplitude (B) and velocity (C) all decreased along the intestine. However, the rates of decrease 
of each were different, as illustrated in the graphs. Quantitative values are included in Appendix. 



'REGROUPS' algorithm for grouping ATs into coherent wave 
fronts were found to require extensive manual review and 
correction. Spatial analysis showed slow wave propagation to be 
dynamic, with propagation occurring in the antegrade and retro- 
grade directions (Figs. 3 A and 3B). Collisions were observed be- 
tween slow waves propagating in opposite directions (Fig. 3C). 
Propagation patterns were dynamic in time, with observed 
changes from antegrade to retrograde propagation within the du- 
ration of the recording periods. Of the 24 total recordings ana- 
lyzed, 8 (33%) exhibited solely antegrade propagation across the 
entire array, 3(13%) exhibited solely retrograde propagation, and 
13 (54%) exhibited both antegrade and retrograde propagation 
(e.g., colliding wavefronts). 

Analysis of Slow Wave Characteristics 

The frequency, amplitude, and velocity of slow waves all de- 
creased along the length of the intestine from proximal duode- 
num to terminal ileum, but the pattern of decrease varied between 
each characteristic, as shown in Figure 6. Slow wave frequency 
gradually decreased in a linear fashion from proximal duodenum 
to terminal ileum. Slow wave amplitude decreased very rapidly 
over the short segment of intestine from the proximal duodenum 
to the proximal jejunum (approximately 0.5 m), and then stabi- 
lized from there to the terminal ileum. Slow wave velocity de- 
creased from proximal duodenum to mid-jejunum, and then re- 
mained relatively stable from there to the terminal ileum. No stat- 
istical differences were observed between the frequency, ampli- 
tude or velocity of slow waves propagating in the antegrade ver- 
sus retrograde direction, so results are reported independent of 



propagation direction. These frequency, amplitude and velocity 
trends are similar to those documented in other species,"^'^'^"^ fur- 
ther validating the analysis pipeline developed in this paper. 

Discussion 

This study presents a platform for obtaining reliable small in- 
testine slow wave recordings using mass-produced flexible arrays 
and computational algorithms for automated slow wave analysis. 
Several new or improved approaches to small intestine mapping 
were developed, notably: (1) the flexible electrode arrays allowed 
for mapping around the intestinal circumference and greater than 
2.5 times the length of intestine mapped in previous studies,^'^^ 
(2) the effect of digital filtering methods were quantified, (3) au- 
tomated algorithms for small intestine slow wave AT marking 
were optimized and (4) the data analysis and visualization process 
was largely automated, including the calculation of slow wave fre- 
quency, velocity and amplitude, and the generation of AT, veloc- 
ity field and amplitude maps. In combination, these outcomes 
constitute an improved platform for small intestine HR mapping 
that will allow more thorough serosal mapping, as well as sub- 
stantial efficiency gains in the data analysis process, while also 
providing enhanced opportunities for clinical application. 

The findings of this study match the trends of decreasing fre- 
quency, amplitude, and velocity along the intestine from histor- 
ical data of other species, particularly the canine (which is the pri- 
mary model for historic in vivo large animal studies),"^'^'^^ provid- 
ing experimental validation for the more comprehensive methods 
developed by this study. However, some quantitative differences 
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exist between the porcine and canine slow waves. The porcine 
small intestine was much longer than that of the canine 
(approximately 15.5 m porcine versus approximately 3 m ca- 
nine"^'^), potentially contributing to some of the interspecies dif- 
ferences in quantitative slow wave results. The decrease in slow 
wave frequency along the porcine intestine was similar to that de- 
scribed in several previous canine studies. ^'^'^^ Although the duo- 
denal frequencies were similar, the porcine terminal ileal fre- 
quency was less than that of the canine (8.4 cpm porcine versus 
12.2 cpm canine^), potentially resulting from the longer porcine 
intestinal length. Porcine duodenal velocity was much lower than 
that of the canine (14.9 mm/sec porcine versus 105 mm/sec ca- 
nine^), although velocity in the terminal ileum was similar (8.8 
mm/sec porcine versus 8 mm/sec canine^). Slow wave amplitude 
is affected by the size and material of the recording electrode, 
making comparisons between studies using different recording 
methods difficult and unreliable. 

Propagation pattern was more variable in the porcine intes- 
tine when compared to a previous HR canine mapping study; ^ 
54% of porcine recordings displayed both antegrade and retro- 
grade propagation (e.g., colliding wavefronts), whereas only 2% 
of the canine recordings displayed such activity, termed 
'propagation blocks' in that study. ^ This increase in wavefront 
collisions could be due to the increased length of the electrode ar- 
ray used in this porcine study (124 versus 46 mm in the canine 
study'), increasing the likelihood that a collision site would be lo- 
cated within the mapped area. However, the increase in multiple 
propagation directions may also be indicative of dysrhythmic slow 
wave activity, which can occur under these experimental con- 
ditions in weaner pigs,^^ further validating the ability of the de- 
scribed methods to map complex propagation patterns. 

The techniques applied in this study can readily be translated 
to human studies, because PCBs can be sterilized and used in- 
tra-operatively.^'^ Unlike in the stomach, where non-invasive cu- 
taneous electrogastrography provides a basic method of evaluat- 
ing gastric slow wave abnormalities,^^ cutaneous intestinal analy- 
sis yields signals of very low signal-to-noise ratio that are difficult 
to interpret and are rarely used.^^ Consequently, there is a critical 
lack of studies on human small intestine slow wave patterns, espe- 
cially relating to the pathology of intestinal motility disorders, 
which remain a considerable clinical challenge.^^ The clinical 
drive for an improved understanding of small intestine slow wave 
events has recently been reinvigorated by a growing under- 
standing of ICC pathologies. For example, the mutation of the 
Nayl .5 sodium channel has recently been associated with irritable 



bowel syndrome. The methods described in this study also cre- 
ate a means to investigate the bioelectrical implications of motility 
disorders like diabetic dysmotility,^ chronic idiopathic pseu- 
do-obstruction,^^ and mitochondrial neurogastrointestinal ence- 
phalomyopathy.^^ 

The recordings obtained and analyzed in this study required 
invasive techniques, as will further investigations using similar 
techniques in the near future. However, the placement of electro- 
des directly on the serosa offers the opportunity to map slow 
waves in the highest detail in normal and diseased human 
intestine.^'^ Once abnormal clinical patterns are defined using a 
HR electrode array as described in this study, a less-invasive de- 
vice with fewer channels could potentially be employed to detect 
dysrhythmia based on signatures of the dysrhythmic patterns, 
such as changes in amplitude and velocity of extracellular 
signals. In the future, once dysrhythmias are better defined by 
methods like those described in this study, much smaller devices 
can be utilized and inserted endoscopically, as have already begun 
to be developed.'^^ While much of the intestine may be un- 
reachable by endoscope, the duodenum is easily accessible endo- 
scopically, and duodenal electrical abnormalities have been im- 
plicated in motility disorders.^^ 

The small intestine is topologically a flexible cylinder held by 
the mesentery. It decreases in circumference along its length and 
is typically collapsed in the fasted state. These attributes, com- 
bined with the biological variability between animals, make it 
challenging to achieve good electrode contact between a HR elec- 
trode array and the intestinal serosa, especially in vivo. The 75 ± 
9% contact achieved in this study is a reflection of the 3-dimen- 
sionally complex organ geometry (effectively a tube that is held 
into curved 'loops' by the mesentery), but this electrode contact 
proved sufficient for spatial analysis of complex and variable slow 
wave propagation patterns. Importantly, the flexibility of the elec- 
trode arrays allowed mapping around the intestinal circum- 
ference for the first time and covered approximately 80% of the 
circumference in the jejunum. 

Appropriate filter settings are vital for the analysis of gastro- 
intestinal slow wave signals, and the comparison of digital filters 
performed in this study is a valuable analysis. The specified 
Savitzky-Golay filter was determined to be the preferred filtering 
method because it maintained the slow wave amplitude and gra- 
dient, which were deemed the most important factors, as only 
mild noise correction is generally necessary and maintaining sig- 
nal quality is vital. The Butterworth 120 cpm filter performed 
similar to the Savitzky-Golay, but the Savitzky-Golay is a time 
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domain based filter that better maintains the slow wave peaks, as 
demonstrated by the amplitude maintenance.^^ 

Intestinal HR mapping yields a vast data set, creating a sub- 
stantial data analysis throughput problem. Slow wave detection, 
analysis, and visualization have previously been performed man- 
ually in the small intestine, ^'^ presenting a time-intensive barrier 
to experimental and translational applications. An automated al- 
gorithm for marking slow wave activity in conscious animals has 
previously been presented by Lammers et al,^^ and is very similar 
to the amplitude-sensitive differentiator approach tested in this 
study. The application and tuning of the FEVT algorithm pre- 
sented in this study constitutes an improvement over this pre- 
vious approach, and the addition of automated frequency, am- 
plitude, velocity and visualization techniques yields an invaluable 
pipeline of analysis tools. 

The FEVT parameters summarized in Table 1 represent the 
set of parameters that yielded the best combined performance 
across all 3 validation data sets, although they were not the opti- 
mal parameter combination for any of the individual data sets. 
Achieving a PPV and Sensitivity of 90% and 87%, respectively 
represents a functional automated detection algorithm that is suf- 
ficiently accurate, although manual review will still be necessary 
when higher reliability is desired. 

Regenerative actions potentials, termed 'spikes,' are also 
known to occur in the small intestine and are closely related to 
contractile response.^ Spikes usually occur during the plateau 
phase of the slow wave, but are also known to occur in the absence 
of slow waves, likely modulated by neural inputs, and resulting in 
propagating contractions. ^'^^ Spikes can propagate independent of 
slow waves, and are characterized by different velocity and prop- 
agation patterns than the corresponding slow wave activity.^^'"^"^ 
Although it was not a focus of the current study, the experimental 
methods developed here could be adapted to map spike activity in 
the future, and a new automated method for spike detection is be- 
ing developed. 

The recordings in this study were performed on fasted, anes- 
thetized pigs, shortly after laparotomy. The effects of anesthesia 
on slow wave activity have not been extensively studied. Howev- 
er, previous mapping studies have been successfully performed 
under anesthesia in several species. ^'^'^^ To mitigate adverse ef- 
fects in this study, intestinal handling was kept to a minimum, 
and temperature and moisture levels were controlled. 

In the future, automating the process of detecting and delet- 
ing the channels that do not contain slow wave signal content, po- 
tentially through frequency or kurtosis analysis, would further 



streamline the data analysis pipeline and eliminate the need for 
these channels to be manually screened. Additionally, automated 
clustering of slow wave events from individual electrode channels 
into corresponding wavefronts requires improvement. Although 
the 'REGROUPS' algorithm was utilized in this study and per- 
formed adequately on selected recordings of simple prop- 
agation,^"^ extensive manual correction was required when prop- 
agation was complex. 

In conclusion, the methods developed in this study allow for 
improved and efficient in vivo studies of small intestine slow wave 
activity. They represent an overall improved pipeline of experi- 
mental recording to data analysis to visualization. Importantly, 
the methods are applicable for use in human studies. Future 
studies will now apply these methods to better define organ-level 
mechanisms of normal and aberrant small intestine slow wave ac- 
tivity, including in human subjects, aiming to improve the clinical 
understanding of small intestine motility in health and disease. 
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Appendix. Average Slow Wave Characteristics Along Intestine 



Location 



Distance from pylorus (m) 



Frequency (cpm) 



Amplitude (|iV) 



Velocity (mm/sec) 



Proximal duodenum 
Proximal j ej unum 
Mid-jejunum 
Proximal ileum 
Terminal ileum 



0.1 

0.5 
4.5 
9.5 
15.5 



16.6 ± 0.4 
16.4 ± 0.3 
13.8 ± 0.5 
11.6 ± 0.8 
8.4 ± 0.7 



145 ± 22 
43 ± 6 
32 ± 6 
29 ± 2 
42 ± 9 



14.9 ± 1.2 
13.4 ± 1.7 
8.8 ± 0.6 
7.6 ± 0.3 
8.8 ± 0.7 



Frequency, amplitude and velocity at 5 discrete locations along the intestine. Frequency, amplitude and velocity all decreased along the length of the intestine, in accord 
with the trends shown in Figure 6. 
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